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ABSTRACT 

We present broad band imaging observations obtained with the "peak up" imagers of the Spitzer 
Space Telescope Infrared Spectrograph (IRS) at wavelengths of 16 fim and 22 fim for a number of 
sources detected primarily at submillimeter wavelengths, which are believed to be at high, though 
undetermined, redshift. We targeted 11 sources originally detected by SCUBA and 5 submillimeter 
sources detected at 1.2mm by MAMBO. Two optically discovered quasars with z>6 were also observed 
to determine if there is detectable dust emission at such high redshifts. Seven of the submillimeter 
sources and both high-rcdshift quasars were detected, and upper limits of about ^50 /iJy apply to the 
remainder. Using their mid-/far-IR colors, we demonstrate that all of the submillimeter sources are at 
z>1.4. The mid-IR colors for two of our detections and three of our strong upper limits suggest that 
these galaxies must be at z>2.5. We also introduce a technique for estimating redshifts and source 
characteristics based only on the ratio of fluxes in the 16 fim. and 22 /im images. 

Subject headings: dust, extinction — infrared: galaxies — galaxies: submillimeter — galaxies: high- 
redshift - galaxies: AGN — galaxies: starburst 



1. INTRODUCTION 

It is currently widely accepted that the energy produc- 
tion in the universe, often quantified with the star forma- 
tion rate (SFR), increases by more than an order of mag- 
nitude as we move from the local Universe to z>2 (e.g., 
iMadau et al.1 11998'). Interstellar dust, a direct product 
of star formation and evolution, often obscures and ab- 
sorbs the light emitted from stars or non-thermal sources 
and reradiates at longer wavelengths in the mid- and far- 
infrared (mid-IR/far-IR) and submillimeter. The results 
of IRAS, launched more than 20 years ago, have demon- 
strated that in the local universe 6% of the total infrared 
emission originates from the centers of galaxies which are 
highly enshrouded by dust and therefore are extremely 
difficult to probe through direct observations in the opti- 
cal and near-infrared (near-IR). These are the luminous 
and ultraluminous infrared galaxies (LIRGs/ULIRGs) 
(e.g.,|Houck ct al. 1984). Deep IR extragalactic surveys, 
mainly performed at 15/zm with the Infrared Space Ob- 
servatory (ISO), recently revealed a strong evolution in 
galaxy counts such that at Z'^l a large fraction of the 
cosmic infrared background (CIRB) is dvie to LIRGs and 
ULIRGs (Franccschini et al. 2001; Elbaz et al. 2002a; 
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lElbaz fc Cesarskvl l2003': 'Me tcalfe et all 1200?). Further- 
more, deep maps with the Sub-mm Common User Bolo- 
metric Array (SCUBA), as well as studies with the Max- 
Plank Millimeter Bolometer (MAMBO), have shown that 
almost 100% of the submillimeter background can be re- 
solved into individual source s, most of which are too faint 
for optical spectroscopy fsee lBlain et al. ir2002l for a re- 
view) . 

Unlike the majority of ISO sources detected at 15/im, 
for which direct optical im aging and redshift mea- 
surements were fe asible (see iFranceschini et aLl 1200 It 
lElbaz et al.l l2002aj) . the nature of the submillimeter 
galaxy (SMG) population is much more elusive. More 
than 300 SMGs have been detected to date, but their 
faintness in the optical and near-IR and the positional 
uncertainty of the '^13" SCUBA beam made a secure 
identification of a counterpart at those wavelengths very 
challenging. The typical spectral energy distribution 
(SED) of a dusty galaxy produces a negative K correction 
in the mm and sub-mm for high redshift systems, mak- 
ing detection of such sources much easier than at shorter 
wavelengths. The observed surface density of SCUBA 
sources and estimated photometric redshifts imply that 
the co-moving density of a high redshift, ~ 10^'^ L© pop- 
ulation of galaxies is ~400 times greater than in the 
local universe (Bl ain et al. Il2003) . Confirmation that 
such high-redsh i ft sou rces are present was achieved by 
iChanman et al.1 l)2003() who used optical spectroscopy 
and accurate radio continuum positions fB arger et alJ 
2000; Ivison et al. 2002) to measure the redshifts for ten 
such sources and show that for this population the me- 
dian redshift is < 2; >~2.4. 

It remains crucial to develop diagnostics for under- 
standing the nature of these luminous, primordial sources 
and to determine accurate redshifts so that their evo- 
lution can be traced. Optical observations probe the 
rest-frame ultraviolet or blue continuum, which only 
characterizes the small fraction of luminosity that es- 



Charniandaris et al. 



capes the dust. Observations at longer rest frame wave- 
lengths can place important constrains on the properties 
of dust as well as the physics of the dominant ionizing 
source and are es sential if the optical is invisible (i.e. 
iFraver et al.ll20n4l . Near-IR spectroscopy even with 8m 
telescopes is extremely difficult for SMGs, which typi- 
cally have an mj > 22mag, so the mid-IR is the next 
available wavelength window. Combining the results 
of the deep surveys performed with ISO and SCUBA 
over the past 7 years, only 11 galaxies have been ob- 
served at both 15^m and 850/im, and redshifts have 
been determined for 9 of them (see Fig. la). Seven 
of the SMGs were detected because of favorable grav- 
itational amplifi cation of their emiss ion by foreground 
lensing clusters (iMetcalfc ct al."2003"), two were identi- 
fied in the CFRS14 field (LiUv et al. 1999: Flores et all 
[1999; Hie don et alJl2004D one was found in the H ubble 
Deep Field IjHuehes et alJll99llAussel. et al.lll999H . and 
one, HRIO, is an ex tremely red object (ERO) at z—lAA 
(|Elbaz et alJ l2002b'). 

The superb sensi tivity of the S pitzer Infrared Spectro- 
graph (IRS) (Hou cket al.ll2004|) . and the ability of the 
IRS to efficiently obtain 16/im and 22/im images to lev- 
els of several tens of fiJy and follow-up spectra to a level 
of ~0.5mJy opens a unique opportunity to examine the 
mid-IR emission from SMGs. In this paper we present 
our first results from a program of targeted observations 
of a number of such sources. We discuss selection criteria 
and observing strategy in section 2, results are presented 
in section 3, and our conclusions can be found in section 
4. 

2. OBSERVATIONS 

The targets were selected from a number of recent sub- 
millimeter catalogs, with several objectives. Primarily, 
we wished to obtain mid-IR images of sources for which 
no positive identification had been made with an optical 
counterpart, or for which no redshift had been obtained 
but was expected to be z>l. Various sources with no op- 
tical counterparts will be found by Spitzer deep surveys 
in the mid-IR, so it is important to begin calibrating the 
properties of samples derived in other ways which may 
have similarities to such sources. The 16/^m to 850/im ra- 
tios can improve estimates of the probable redshift. We 
can also utilize the 16/im to 22/im flux ratio to deter- 
mine the slope of the infrared spectrum or the presence 
in one or the other band of any strong polycyclic aro- 
matic hydrocarbon (PAH) emission features or silicate 
absorption features. The final objective of our imaging 
survey of SMGs was to determine whether any are bright 
enough for an IRS spectroscopic follow-up which can give 
a precise redshift measurement^. In addition to observ- 
ing the submillimeter sources, we added two quasars at 
very high redshift (z > 6) identified by the Sloan Digi- 
tal Sky Survey to determine if such distant quasars also 
have sufficient re-emission from dust to be detectable in 
Spitzer deep surveys. One of these quasars, J1148-I-5251, 
has many characteristics of the extreme submillimeter 
sources; it has ~ 2 x 10^ ° Mp, of molecular gas and an 
L(FIR)- 1.3 X 1013 Lq llWalter et alJ 12001 . has been 
detected by MAMBO (fi.2mm=5.0±0.6mJy), and is es- 

^ For f22fim >0.7mjy wc can obtain a 20-38/ini spectrum with 
a > 5(T detection in ~lhr. 



timated to contain ^ 7 x 10^ M© of dust (|Bertoldi et alJ 
12003(1 . Whether sufficiently hot dust is present in such 
AGN to be seen at the much shorter wavelengths ob- 
served in the rest frame is important to know for inter- 
preting other sources. All of our targets were observed 
during the first 3 IRS campaigns that took place between 
14 December 2003 and 9 February 2004. 

We used the blue and red peak up cameras located in 
the short low module of the IRS to observe our targets 
at IGfim (13.5-18.7/im) and 22 /im(18.5-26.0/im). By 
using the Spitzer focal plane as the reference frame, it 
is possible to program observations to obtain dithered 
observations in both the 16/Ltm and 22 /Ltm images as a 
stan dard IRS short lo w first order (SLl) staring AOT 
(see lHouck et al.|l2004() with the targets defined of type 
"fixed cluster offsets" on the plane of the array. The off- 
sets, in arcsecs parallel and perpendicular to the SLl slit, 
were set to (-62", 289") and (50", 284") for the centers 
of the blue and red peak up windows respectively. When 
the telescope was performing the standard nodding pat- 
tern on the SLl slit, the actual target was imaged at two 
positions on the peak up cameras separated by 18". This 
approach enables us to obtain both images with a single 
telescope pointing, and also to collect the data in the 
"raw mode", which allows flexible image analysis. 

The data were processed by the standard IRS pipeline 
(version S9.1) at the Spitzer Science Center (see chap- 
ter 7 of Spitzer Observing Manual^). The 2D images 
were converted from slopes after linearization correction, 
subtraction of darks, and cosmic ray removal. The re- 
sulting images were divided by the photometric flat, and 
a world coordinate system was inserted into them using 
the reconstructed pointing of the telescope. The astro- 
metric accuracy of our images is better than ^2" and 
the FWHM of the point spread function is 3.5" and 4.8" 
at 16/im and 22/im respectively. Since we typically had 
multiple images of the same target, we used the astrom- 
etry information to shift and co-add our pointings, as 
well as to perform sky subtraction by taking the differ- 
ences of the two nod positions. Flux calibration was per- 
formed using a number of stars for which high quality 
spectral templates were available. The properties of all 
targets observed are presented in Table 1. Our photomet- 
ric uncertainties are less than '^10%. The detection limit 
of the peak up cameras depends rather strongly on the 
brightness of the background. Since the targets observed 
are in a low background region with surface brightness 
^20MJysr^i, the la detection limit is ~30/iJy for an 
integration time of 120sec at 16/im and 240sec at 22/im. 

3. RESULTS 

As we show in Table 1, we were able to detect 10 out 
of the 16 submillimeter targets at > 3a, including all 3 
sources for which there are published r edshifts. Only on e 
source, CFRS 10.1411 at just z=0.07 fLi llv et alJ ll999'l. 
is above ImJy. The fluxes of the remainder of the SMG 
detections are on the order of a few hundred /i Jy in both 
16 and 22/im. No source without a known redshift is 
above the f22Ami=0.7mJy. The two quasars, J1048-f 4637 
and J11484-5251, with redshifts of 6.23 and 6.42, respec- 
tively, were also detected. Given their extreme redshift, 
our measurements probe their rest frame near-IR emis- 

^ http://ssc.spitzcr.caltcch.edu/docunicnts/som/ 
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sion ('^2-3/ini) which is hkely due to a hot dust com- 
ponent. This is especially significant for J1148-I-5251, a 
strong submillimeter source, because it could be used to 
explore whether AGN at high redshift which have suffi- 
cient cool dust for submillimeter emission can contain a 
hot dust component. It also illustrates the challenges in 
determining the primary luminosity source (starburst or 
AGN) for the "submillimeter" population. 

When wc combine our measurements from the IRS 
peak up filters with the SCUBA 850/ini points, we can 
bracket the peak of the infrared SED of a source, typi- 
cally found at ~80/xm in the rest frame, and obtain a pho- 
tometric estimat e of its redshift. I n a manner similar to 
the discussion of iLillv et alJ l|1999() in exploring the ISO 
15/im to 850/im colors, we calculated the ii6fj.m/is50fj.m 
ratio as a function of z for a number of various tem- 
plate SEDs. Our results for a few selected templates are 
presented in Figure la. We can see that as the source 
moves to higher z, due to the negative K correction the 
flux sampled by SCUBA at 850/im remains relatively 
flat or decreases only slightly. The rate of the decrease 
clearly depends on where the far-IR SED peaks and it is 
faster (slower) when the galaxy has warmer (colder) far-IR 
color. The flux measured by our 16/Ltm peak up window 
decreases much more rapidly up until z~l when the rest 
frame ~6-8/im range begins to enter the window. Be- 
yond that point the details depend on features in the 
mid-IR spectrum of the source. If a source has a feature- 
less mid-IR continuum with weak PAH emission such as 
an AGN dominated template (i.e. NGC 1068), the fig/im 
flux would decrease monotonically. The situation is a bit 
more complex if strong PAH emission features (such as 
those seen in M82) or silicate absorption (as in the case of 
Arp 220) are present in the spectrum. It is evident from 
Figure la that given a sufficiently large sample of possi- 
ble templates we can define a lower limit to the possible 
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at a given redshift. Calculating this ratio as 
a function of z for a wide range of known SEDs suggests 
that at any given z>0.5, where most SMGs are found, an 
SED similar to that of Arp 220 provides the most con- 
servative limit of the possible value fie/jm/fsso/im- By 
measuring this ratio for our SMG targets, we can obtain 
the absolute lower limit for their redshift, assuming this 
SED. 

To demonstrate this, we also plot in Figure la, along 
with our data, all SMG sources for which redshifts and 
ISO 15/im fluxes have been published to date^. We see 
that all SMGs of known redshift have fie/im/fsso/jm ra- 
tios placing them between Arp 220 and Mrk 273. An ac- 
curate estimate of the redshift clearly depends on the 
assumed SED of the source. Our conservative method 
indicates that with one exception^ all SMGs in Table 
1 with no known redshift must be at z>1.4. Even if at 
Z'^l.S, a source like Arp 220 would have an mj ~23.5mag, 
well within spectroscopic capabilities from the ground. 
Our upper limits for the fie^m/fssopm ratio on sources 

® Given the similarities between the ISO IS^m (LW3) and the 
IRS 16/ini blue peak up filters, as well as the possible shapes of the 
mid-IR SEDs for the ISO detected SMGs we have set their fi6/jm 
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^ We cannot explain at this point why source #9 despite its 
strong det ection and the fact that radio/sub-mm data suggest a 
z~2.4 iBa rger et aLIEOOfl^ . has rather peculiar mid-IR and sub- 
mm colors. 



#4, #10 and #14 suggest that these sources are beyond 
z>2.5. 

Interestingly, the highest redshift (z=2.8) know n ISO 
sourc e, SMM 02399-0136 in AbeU 370 (Metcalf e etal J 
31). as wefl as our target #1, the ERO N4 of 
Fraver et al.l ll2003t ). are located very close to the pure 



AGN locus (NGC 1068). If the SMGs are indeed warm 
AGN-like sources following the SEDs discussed bv lBlainI 
(1999), then our Arp 220 limit is extremely conservative. 
In this case, the luminosity of these sources would derive 
from a scaled up version of an enshrouded AGN, and the 
thermal emission from the AGN in the 1-6/im rest frame 
will enter into our 13/im-18/iin blue peak up window at 
z>2. This will increase the expected fie^m/fssopm ra- 
tio by a factor of ~8, as seen in Figure la, so that the 
necessary redshifts will easily exceed z~ 3. 

We can also begin to investigate whether it will eventu- 
ally prove possible to get reliable photometric redshifts 
from the 16/im and 22/im imaging data alone, and to 
use this mid-infrared flux ratio for diagnostic informa- 
tion on the dominance of an AGN or a starburst. The 
concept of this is illustrated in Figure lb. Galaxies 
whose mid-IR luminosities derive primarily from star- 
bursts are characterized by s trong PAH emission fea- 
tures (e.g. iBrandl et al.l 12013^ . whereas AGN generally 
produce dust continua with out PAH emission b ut often 
with sihcate absorption fe.g. lLaurent et aLll200nD . While 
we intend to calibrate the effect with many more spec- 
troscopic observations of known starbursts and AGN, we 
illustrate in Figure lb the changes that would be ex- 
pected in the fi6/im/f22/jm ratio for templates like the 
ISO spectra of the starburst M 82 (Stu rm et al.l2000(l . the 
AGN NGC 1068 l|Le Floc'h et al...200l. a nd the more 
enshroude d sources Arp220 ((SDOon et a l. 2004) and 
Mrk 273 ( Devriendt et al.lll999fl . Note that even though 
NGC 1068 hosts a Sy2 nucleus affected by considerable 
extinction in the optical, it displays a nearly power-law 
continuum w ith very weak PAH and silicate features i n 
the mid-IR (|Sturm et all l2000l iLe Floc'h etUI l200l|) . 
Over a broad range of redshifts for z>l, the approxi- 
mately power-law spectrum of this AGN means, there- 
fore, that the 16/im to 22//m flux ratio has a narrow 
range, from 0.7 to 0.9, whereas this ratio for other sources 
with stronger spectral features ranges between extremes 
of 0.1 to 1.7. It is clear that this flux ratio information 
is potentially useful for constraining source characteris- 
tics and redshifts. For illustration, we apply this con- 
cept to the present results in Table 1 by plotting sources 
in Figure lb at the redshifts corresponding to the min- 
imum redshift determined from the iiep,m/i850p,m ratio. 
Only sources #6, #9, and #11 have ratios considerably 
outside the expected NGC 1068 values. Having the ra- 
tios shown implies that these 3 sources have strong spec- 
tral features and, for sources #6 and #11, yields red- 
shift estimates consistent with those derived from the 
fi6;xm/f850;xm ratio in Figure la. 

4. CONCLUSIONS 

The unprecedented sensitivity of the IRS peak up cam- 
eras has enabled us to perform efficient imaging of 16 
SMGs at 16/im and 22/im. Ten of these sources are 
detected, which doubles the number of mid-IR detec- 
tions for sources discovered through submillimeter sur- 
veys. Using a number of SED templates along with the 
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TABLE 1 
MiD-IR Properties of high-z sub-mm Galaxies 



RA DEC 

(J2000) (J2000) 



/i6m 



jKi^m J22fj.m t 

(mjy) (mjy) 



/850/jm 

(mJy) 



fl6iim 



[ID] /Comments 



04:43:07.1 


02:10:25.1 


0.52 


0.98 


0.54 


2.51 


7.2 


0.072 


1 


SMM J0443H-0210(N4) in (1) 




10:00:37.2 


25:14:57.9 


1.88 


2.55 


0.73 


0.07 


4.7 


0.400 


2 "" CFRS 10.1411 in (2) 




10:00:38.2 


25:14:49.2 


0.27 


0.62 


0.42 


0.50 


4.8 


0.056 


3 CFRS 10.1153 in (2) 




12:36:06.4 


62:11:38.0 


<0.026 


<0.030 


- 


>2.5 


15.4 


<0.002 


4 HDFSMM-3606+1138 in (4) 




12:36:18.4 


62:15:49.2 


0.12 


0.13 


0.92 


>1.5 


7.8 


0.015 


5 #7 in (3) 




12:36:46.7 


62:14:46.9 


0.13 


0.38 


0.34 


>1.5 


10.7 


0.012 


6 #33 in (3) 




12:36:52.32 


62:12:26.3 


<0.025 


<0.041 


- 


>2.0 


7.0 


<0.004 


7 HDF N(850.1) in (5) 




12:37:00.4 


62:09:08.5 


0.07 


0.13 


0.54 


>1.8 


11.9 


0.006 


8 #49 in (3) 




12:37:09.5 


62:08:39.5 


0.70 


0.61 


1.15 


>0.5 


8.3 


0.084 


9 #58 in (3) 




12:37:29.7 


62:10:51.0 


<0.028 


<0.024 


- 


>2.5 


14.3 


<0.002 


10 


HDFSMM-3730-f 1051 in (4) 




14:17:40.2 


52:29:06.5 


0.19 


0.14 


1.35 


>1.4 


8.8 


0.022 


11 


CFRS 14A in (6) 




12:05:17.59 


-07:43:11.5 


<0.060 


<0.068 


- 


>1.6 


6.3 


<0.009 


12 


" MMJ 120517-0743.1 (NTT-MM25) in 


(7) 


12:05:39.45 


-07:45:26.1 


<0.049 


<0.050 


- 


>1.7 


6.3 


<0.008 


13 


^ MMJ 120539-0745.4 {NTT-MM16) in 


(7) 


12:05:46.55 


-07:41:33.2 


<0.040 


<0.060 


- 


>2.2 


18.5 


<0.003 


14 


'^ MMJ 120546-0741.5 {NTT-MM31) in 


(7) 


15:41:25.2 


66:14:39.6 


0.11 


0.16 


0.69 


>1.6 


10.7 


0.010 


15 


J 1541-1-6615 {A2125-MM26) in (7) 




15:41:27.3 


66:16:59.4 


0.05 


0.07 


0.71 


>2.2 


14.6 


0.003 


16 


J 1541-1-6616 (A2125-MM27) in (7) 




10:48:44.9 


46:37:15.5 


0.23 


0.37 


0.62 


6.23 






17 


Jl048-h4637 




11:48:16.5 


52:51:49.1 


0.51 


0.74 


0.69 


6.42 






18 


J1148-I-5251 





Note. — In cases of non detections we ind icate our la upper limit s. T he identificat i ons an d m ore detail on in dividual sourc es are presente d in: 
flilFraver et al.l 1200311 . (•2') ILillv et"!!! 1199911 . O'l lBarger et al.l i200Cll . f4) IBorvs eTaT] i2002f) . fS'l lHughes et al.l 119981. (6) E alcs et al.H1999ll . (7) 
lEales et al.l 1200311 . Redshifts in brackets arc absolute minimum redshifts based on photo-z and an Arp 220 template (see section 3 and Fig.l). 

^Unlike the rest of the sources where the on source integration times were 120sec and 240sec at 16;im and 22/im, the corresponding times for these 
targets were 140scc and 280scc respectively. 



observed ii6fj.m/ia50fim ratios or upper limits, we have 
shown that 13 of these SMG sources must be located 
at z>1.4, with some being constrained to z>2.5. We 
also introduce use of the observed fi6^m/f22/im ratio to 
constrain photometric redshifts, and to demonstrate that 
some of the SMGs must have mid-infrared spectra similar 
to those of AGN. 
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Infrared Imaging of High Redshift Submillinictcr Sources 
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Fig. 1. — a) The expected ratio of the flux densities fie/jm/fsso/jm of a distant source as a function of redshift if the SED of the source 
is similar to a starburst Uke M82, an AGN Hke NGC 1068, or an enshrouded source such as Arp220 or Mrk273. Note that the SED of 
Arp 220 defines in the most conservative way the minimum redshift of the SMG source for a given value of the ratio. The minimum redshift 
for each SMG of our sample is indicated by the source number from Table 1. All previously known ISO 15/im sources at z>l for which 
both redshifts and 850/im fluxes are available are plotted as asterisks. Interestingly, the most distant ISO source as well as our source 1 
for which the redshift is known, fall most closely to the AGN template. Six sources (4,7,10,12,13, and 14) for which we only have upper 
limits at 16 and 22fim are indicated with arrows and must be at even higher redshifts (sources 4,10, and 14 must be at z>2.5). b) The 
expected flux ratio fi6/jm/f22/jm as would be measured with the IRS peak up filters for the templates used above. Note that if a galaxy 
has a starburst SED, similar to M82 which is dominated by PAH emission in the 6-8/im rest frame, or strong silicate absorption, such 
as Arp 220, then the fi6/jm/f22/jm ratio may vary up to a factor of 8 as these features pass in and out of the filters, and this ratio can 
potentially be used as a tracer of the redshift of the source. The measured ratios for our sources plotted against their minimum possible 
redshift are indicated by the source ID (see Section 3). 



